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The recent past witnesses the growing interdependence between the physics of hadrons, the 
physics of relativistic heavy-ion collisions, and the physics of compact objects in astrophysics. A 
notable example is the kaon which plays special roles in all the three fields. In this paper, we 
hrst review the various theoretical investigations of kaon properties in nuclear medium, focusing on 
possible uncertainties in each model. We then present a detailed transport model study of kaon 
production in heavy-ion collisions at SIS energies. We shall discuss especially the elementary kaon 
and antikaon production cross sections in hadron-hadron interactions, that represent one of the most 
serious uncertainties in the transport model study of particle production in heavy-ion collisions. The 
main purpose of such a study is to constrain kaon in-medium properties from the heavy-ion data. 
This can provide useful guidances for the development of theoretical models of the kaon in medium. 

In the last part of the paper, we apply the kaon in-medium properties extracted from heavy-ion data 
to the study of neutron star properties. Based on a conventional equation of state of nuclear matter 
that can be considered as one of the best constrained by available experimental data on Hnite nuclei, 
we find that the maximum mass of neutron stars is about 2Mq, which is reduced to about l.SM© 
once kaon condensation as constrained by heavy-ion data is introduced, 
pacs: 25.75.Dw, 97.60.Jd, 26.60.-fc, 24.10.Lx 
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I. INTRODUCTION 


There is currently growing interplay between physics of hadrons (especially the properties of hadrons in dense matter 
which might reflect spontaneous chiral symmetry breaking and its restoration), the physics of relativistic heavy-ion 
collisions (from which one might extract hadron properties in dense matter), and the physics of compact objects in 
astrophysics (which needs as inputs the information gained from the first two fields). A notable example is the kaon 
{K and K), which, being a Goldstone boson with strangeness, plays a special role in all the three helds mentioned. 

Ever since the pioneering work of Kaplan and Nelson ||^,^ on the possibility of kaon condensation in nuclear matter, 
a huge amount of theoretical effort has been devoted to the study of kaon properties in dense matter, based mostly 
on the SU(3) chiral perturbation theory M- Kaons, as Goldstone bosons with strangeness, play a special role in 
the development of hadron models. The mass of the strange quark is about 150 MeV, which is, on the one hand, 
considerably larger than the mass of light (up and down) quarks (~ 5 MeV), but on the other hand, much smaller 
than that of the charmed quark (^ 1.5 GeV). In the limit of vanishing quark mass the chiral symmetry is good, and 
systematic studies can be carried out using chiral perturbation theory for hadrons made of light quarks. Since its 
mass is not negligibly small as compared with the typical QCD chiral symmetry breaking scale A^^ ~ 1 GeV, the 
expansion in terms of the strange quark mass is much more delicate. In spite of this, chiral perturbation calculations 
have been extensively and quite successfully carried out in the recent past for the study of kaon-nucleon (KN) and 
antikaon-nucleon {KN) scattering. 

Part of the reason for this success, at least in free space, can be attributed to the availability of a large body of 
experimental data which can be used to constrain the various parameters in the chiral Lagrangian. In studying kaon 
properties in nuclear matter, a new scale, namely the nuclear Fermi momentum kp, is introduced. This renders the 
chiral perturbation calculation in matter much more subtle than that in free space. It is very important to have 
experimental data available so that the predictions of chiral perturbation calculations can be tested. Kaonic atom 
data provide information on this, but is restricted to very low densities |17 . For the study of kaon condensation, 
densities much higher than that accessible by kaonic atoms are involved. This can only be obtained by analysing 
heavy-ion collision data on kaon spectra and flow. 

Measurements of kaon spectra and flow have been systematically carried out in heavy-ion collisions at SIS (1- 
2 AGeV), AGS (10 AGeV), and SPS (200 AGeV) energies Q. By comparing transport model predictions with 
experimental data, one can learn not only the global reaction dynamics, but more importantly, the kaon properties 
in dense matter. Of special interest is kaon production in heavy-ion collisions at SIS energies, as it has been shown 
that particle production at subthreshold energies is sensitive to its properties in dense matter [pj|-[^. Recently, high 
quality data concerning and K~ production in heavy-ion collisions at SIS energies have been published by the 
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KaoS collaboration at GSI The KaoS data show that the K~ yield at 1.8 AGeV agrees roughly with the yield 
at 1.0 AGeV. This is a nontrivial observation. These beam energies were purposely chosen, such that the Q-values 
(the difference between the available center-of-mass energy and the threshold) for NN NKA and NN NNKK 
are identical and both are about —230 MeV. Near their respective production thresholds, the cross section for the K~ 
production in proton-proton interaction is one to two orders of magnitude smaller than that for AT'*' production . 
In addition, antikaons are strongly absorbed in heavy-ion collisions, which should further reduce the K~ yield. The 
KaoS results of K~ /~ 1 indicate thus the importance of kaon medium effects which act oppositely on the 
and K~ production in nuclear medium, and can thus provide useful information on kaon properties in dense matter. 

Studies of neutron star properties also have a long history, both observationally and theoretically. A recent compi¬ 
lation by Thorsett quoted by Brown shows that well-measured neutron star masses are all less than I.SMq. On 
the other hand, most of the theoretical calculations based on conventional nuclear equation of state (EOS) predict a 
maximum neutron state mass above 2Mq ||2^-|^. The EOS can, therefore, be substantially softened without running 
into contradiction with observation. Various scenarios have been proposed that can lead to a soft EOS. Brown and 
collaborators suggested that kaon condensation might happen at a critical density of 2-4po @|- The existence of 
many K~ mesons implies the coexistence of a large number of protons that are needed to neutralize the negative 
charge of K~ mesons Based on this. Brown and Bethe also proposed the interesting possibility of the existence 
of a large number of low mass black holes in the galaxy . 

In the first part of this paper. Section 2, we review various theoretical approaches for studying kaon properties 
in nuclear medium. We shall focus on the limitations and possible uncertainties in each approach. In the second 
part. Section 3, we present a detailed analysis of kaon production in heavy-ion collisions at SIS energies. We will 
discuss in detail the elementary kaon and antikaon production cross sections in hadron-hadron interactions, which 
are the major sources of uncertainties in the transport model study of particle production in heavy-ion collisions. By 
comparing transport model predictions with experimental data, we try to extract kaon in-medium properties at high 
densities. This information is very useful for the development of chiral perturbation theory for finite density. This also 
provides guidance for the study of kaon condensation in neutron stars, which is the focus of Section 4. We calculate 
neutron star masses based on the nuclear equation of state (EOS) of Furnstahl, Tang and Serot |^l[. Without kaon 
condensation, the maximum neutron star mass is predicted to be about 2.0Mq, which reduces to about IAMq once 
kaon condensation is introduced. The paper ends with summary and outlook in Section 5. 


II. KAON IN DENSE MATTER: A REVIEW 

Ever since the pioneering work of Kaplan and Nelson Q on the possibility of kaon condensation in dense nuclear 
matter, many works have been devoted to the study of kaon properties in nuclear matter. There are two typical 
approaches to this problem. One was initiated by Kaplan and Nelson, and is based on the chiral perturbation theory. 
The other has been pursued by Schaffer et al. and Knorren et al. based on the extension of the Walecka-type mean field 
model from SU(2) to SU(3). In other words, it follows the traditional meson-exchange idea for hadron interactions. 
This latter approach has also been developed by the Jiilich group for the description of kaon-nucleon scattering in 
free space [|32| , p3| . In addition, kaon properties in dense matter have been studied using the Nambu—Jona-Lasinio 
model , treating the kaon as quark-antiquark excitation, and phenomenological off-shell meson-nucleon interaction 
iP. Although quantitatively, the results from these different models are not completely identical, qualitatively, a 
consistent picture, namely in nuclear matter the kaon feels a weak repulsive potential and the antikaon feels a strong 
attractive potential, has emerged. In this section we will discuss mainly the results of the chiral perturbation theory 
and the meson-exchange model. 


A. Chiral perturbation theory 

The interactions between pseudoscalar mesons (pion, kaon, and eta meson) and baryons (nucleon and hyperon) are 
described by the 811(3)^x811(3)^; nonlinear chiral Lagrangian which can be written as 

C = i/^TrS^Sa^Et + i/2A[TrM,(E - 1) + h.c.] + - mB)B 

+ mBY[Vf,, B] + DTiBYl^A^, B} + ETrSy'^ySjA^, B] 

+ aiTri3(^Mq^ -|- h.c.)B + a 2 "TiBB{^Mq^ + h.c.) 

-h a 3 [TrMqE -h h.c.jTiBB. (1) 
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In the above, B is the baryon octet with a degenerate mass ms, and 

S = exp(2i7r//), and ^ =-s/S = exp(i7r//), (2) 

with TT being the pseudoscalar meson octet. The pseudoscalar meson decay constants are equal in the SU(3)v limit 
and are denoted by / = /tt — 93 MeV. The meson vector and axial vector currents are defined as 

and Af, = (3) 

respectively. The current quark mass matrix is given by Mq = diagjmq, mg}-, where we neglect the small difference 
between the up and down quark masses. 

Expanding S to order of 1//^ and keeping explicitly only the kaon field, the first two terms in Eq. ( 0 ) can be 
written as 


d^KdfiK — A{mq + ms)KK + • 


where 


K = 


(S) 


and K = {R- K°), 


and the ellipsis denotes terms containing other mesons. 

Keeping explicitly only the nucleon and kaon, the third and fourth terms in Eq. (|^) become 


3z 


N(i^^d^ — mB)N — dt K + ■ 


where 


N = 


and N = (p n), 


and the ellipsis denotes terms involving other baryons and mesons. 

The last three terms in Eq. (|^) can be similarly worked out, and the results are 


TrBi^Mq^ + h.c.)B = 2mqNN - — (m, + + 


TvBB{^Mq^ + h.c.) = 2msNN -+ ms)KK + 


2NN, 


NN, 

w 

NN, 

7 ^ 


[TrMgS + h.c.jTrSi? = 2{2mq + ms)NN -+ ms)KK + ■ 

Combining above expressions, one arrives at the following Lagrangian, 

C = iV(f7^(9^ - mN)N + - (m^ - ^^NN)KK 


P 


- ^NfNK *dtK A 


where the kaon mass is given by 

m\ = K{mq + ms), 

and the nucleon mass by 

mN = ms — 2[aimq + a2ms + a^{2mq + ms)]. 
Also, the KN sigma term can be expressed as 


(4) 

(5) 


( 6 ) 

(7) 


( 8 ) 


(9) 


( 10 ) 


( 11 ) 
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( 12 ) 


Skat = + ms){N\uu + ss\N) 

= + "*s)(«i + 2 a 2 + 403). 

We note that the last term in Eq. i) is the usual Weinberg-Tomozawa term; it gives rise to a repulsive vector 
potential for K~^, and an attractive vector potential for K~. The one that goes with the kaon mass is usually called 
Kaplan-Nelson term; it provides scalar attraction for both K'^ and K~, and thus reduces their mass. The amount 
of scalar attraction is linearly proportional to the kaon-nucleon sigma term, Tikn- From Eq. ( |l^ we see that Y^kn 
depends on three coefficients, aims, a 2 ms, and a^mg- While the first two coefficients are relatively well determined 
from the baryon mass splitting Q, the dominant term, a^mg, is only poorly known, mainly because of the large 
uncertainties in the nucleon strangeness content. 

On the other hand, the kaon-nucleon sigma term can be related to pion-nucleon sigma term, S,rW @| • The latter 
is relatively well determined to be about 45 MeV p^. 


(jnq + ms){N\uu + ss\N) 
2mq{N\uu-\-dd\N) 

~ t(1 H-)(1 + y)^TrN, 

4 mq 


(13) 


where the strangeness content of nucleon is defined as 


2{N\ss\N) 
{N\uu + dd\N) 


(14) 


From the particle data book we know that the average mass of the light quark is mq « 7.5 MeV. There is a large 
uncertainty in the strange quark mass mg which ranges from 100 to 300 MeV There is also a large uncertainty 

in the nucleon strangeness content, which ranges from 0 to about 0.3 

In the mean-field approximation and including only the Weinberg-Tomozawa and the Kaplan-Nelson terms, the 
kaon dispersion relation in nuclear matter is then given by 


w^(p,Pb) = -bp^ - 


-‘KN 


J2 ' 4 


3 U) 


PN, 


(15) 


where p is the three-momentum of the kaon. From the dispersion relation, the kaon energy in medium can be obtained, 
i.e., 


UK 


-b 


Ykn 

~yrPs + 



3 pn 
87^’ 


and similarly for the antikaon. 


UK = 


‘■K 


Ykn 


PS + 



3 Pn 

sY' 


(16) 


(17) 


Note that the scalar attraction depends on nucleon scalar density pg, which is model dependent |^. 

There are a number of corrections to these simple mean-field results. Below we discuss some of them. 

Range term: At the same order in chiral perturbation theory as the Kaplan-Nelson term is the energy-dependent 
range term. This has been studied in great detail in Refs. [|^,Q. It has been shown that the range term reduces the 
scalar attraction and can be ap pro ximately included, with Ykn ~ 400 MeV, by multiplying the Kaplan-Nelson term 
by a factor 1 — 0.37aj^^/m^ [^. Note that in the nuclear medium, the kaon effective mass increases, while that 
of the antikaon decreases. Therefore, while the range term is very important for K^, it is less so for K~. In other 
words, the scalar attraction for is different from that for K~. 

Brown-Rho scaling: In the nuclear medium, the pion decay constant may decrease js^] . At finite density, this can 
be easily seen from the Gell-Mann—Oakes—Renner relation. In Ref. [Q, Brown and Rho has shown that the scaling 
in /tt is instructive in linking the chiral Lagrangian in free space to Walecka mean-field model for finite nuclei. The 
enhancement of low-mass dileptons in heavy-ion collisions might be considered as empirical evidence for this scaling 
[p9|-pl|. Including this scaling, the vector attraction for K~ increases significantly. 
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Short-range correlations: At high densities, a simple mean-field type treatment of kaon-nucleon and nucleon- 
nucleon interactions might be inadequate. The effects of short-range correlations in nucleon-nucleon and kaon-nucleon 
interactions were studied by Pandharipande, Pethick, and Thorsson ||4^ , and were found to reduce the scalar attraction 
significantly. There will, however, be a tendency for the effects from scaling in and short range correlation to cancel 
each other. 

A(1405) and coupled-channel effects: The isospin averaged K~N scattering length is negative in free space Q, 
implying a repulsive K~ optical potential in the simple impulse approximation. However, a systematic analysis of 
the kaonic atom data shows that the K~ optical potential is deeply attractive, with a value of about 200± 20 MeV 
at normal nuclear matter density Unlike the kaon which interacts with nucleons relatively weakly so that the 
impulse approximation is reasonable, the antikaon interacts strongly with nucleons so we do not expect the impulse 
approximation to be reliable. The antikaon-nucleon (KN) interaction at low-energy is strongly affected by the A(1405) 
which is a quasi bound state of an antikaon and a nucleon in the isospin 1 = 0 channel and can decay into Ett channel 
I p^ . Thus, in principle one needs to carry out a coupled-channel calculation for KN, Att and Stt including the 
effects of A(1405) in both free space and in nuclear medium [p^ , ^4p^ ]. This has been studied in detail by Weise 
and collaborators 115 |Jh| . In effective chiral Lagrangian, however, one can also introduce A(1405) as a fundamental 
field, which is consistent with phenomenological KN scattering amplitudes and branching ratios Because of 

Pauli blocking effects on the intermediate states, the possibility of forming a bound KN state (A(1405)) decreases 
with increasing density, leading to a dissociation of A(1405) in nuclear medium. This induces a transition of the K~ 
potential from repulsion at very low densities to attraction at higher densities j^,^. It was shown in Refs. 
that, because of the diminishing role of A(1405) in dense matter, the prediction for the kaon condensation threshold 
is relatively robust with respect to different treatments of A(1405). 

In view of large uncertainties in and difficulties in treating systematically high-order corrections, we adopt 

in this paper a more phenomenological approach. We assume that the effects from the scaling in and short- 
range correlations approximately cancel each other. Furthermore, we introduce two free parameters, ax and 
which determine the scalar attractions for K'^ and K~. We assume that these are density independent, although in 
principle they should be density dependent, since the range term, the scaling in /,r, and the short-range corrections 
are all density dependent. The kaon and antikaon energy in the nuclear medium can then be written as 


LOk = 


UJx = 


-\-k^ — axPs + {bxPx)^] ^ +hxPN, 

(18) 

1/2 

— axPs + [bxpxT] —bxPN, 

(19) 


where bx = 3/(8/^) « 0.333 GeVfm^. We try to determine ax and ax from the experimental observables in heavy-ion 
collisions. 

Since the kaon-nucleon interaction is relatively weak as compared to other hadron-nucleon interactions, impulse 
approximation is considered to be reasonable for the kaon potential in nuclear matter, at least at low densities. This 
provides some constraint on ax- Using ax = 0.22 GeV^fm^, we find that at normal nuclear matter density po = 0.16 
fm“^, the K'^ feels a repulsive potential of about 20 MeV. This is in rough agreement with what is expected from 
the impulse approximation using the KN scattering length in free space. Note that as p approaches zero. 


ax 


0.63S^jv 

f2 
J TT 


( 20 ) 


Therefore ax = 0.22 GeV^fm^ implies YixN ~ 400 MeV. This is exactly the value determined in Ref. by fitting 
the KN scattering length. 

Determination of the ax is more delicate, as impulse approximation does not apply. The ax should show stronger 
density dependence than ax, as the effects of A(1450) and coupled-channels are both strongly density dependent. We 
will try to determine this value by fitting to heavy-ion data. We find that, when using the model of Furnstahl, Tang, 
and Serot for dense matter, ax = 0.45 GeV^fm^ provides a good fit to the K~ data in heavy-ion collisions at SIS 
energies. Since the K~ production chiefly proceeds at the high densities (2po-3po)j the value of ax determined here 
needs not apply to lower densities, e.g., those sampled in kaonic atoms. 

With these two parameters we show in Fig. the effective masses of kaon and antikaon defined as their energies 
at zero momentum. It is seen that the kaon mass increases slightly with density, resulting from near cancellation of 
the attractive scalar and repulsive vector potential. The mass of the antikaon drops substantially. At normal nuclear 
matter density, the kaon mass increases about 4%, in rough agreement with the prediction of impulse approximation 
based on KN scattering length. The antikaon mass drops by about 22%, which is somewhat smaller than what has 
been inferred from the kaonic atom data |^, namely, an attractive K~ potential of 200 ± 20 MeV at po- 
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From their in-medium dispersion relations, we can define the kaon and antikaon potential as the difference between 
its energies in the medium and in free space , 


Uk,k — ^K.K ~ ( 21 ) 

with ojQ = y/Tn^K + p^- Kaon and antikaon potentials at two different momenta are shown in Fig. The open circle 
in the figure is the K~^ potential expected from the impulse approximation using the kaon-nucleon scattering length 
in free space . The solid circle is the K~ potential extracted from the kaonic atom data . 


B. Meson-exchange model 

The chiral Lagrangian as given above does not describe properly the nuclear matter properties. Therefore in most 
of the studies about kaon condensation in density matter, the nuclear matter properties are obtained based on some 
form of Walecka-type mean-field model; e.g., in this work we use the model of Ref. |^l|, which has imposed chiral 
constraints on the mean field. The question then is, as put forward in Ref. whether it is possible to describe 
the kaon-nucleon interaction using the idea of meson-exchange. This has been addressed by Schaffer et al. for 
symmetric nuclear matter, and by Knorren et al. for more general case as encountered in neutron stars. 

In the meson-exchange picture, the scalar and vector interaction between kaon and nucleon are mediated by the 
exchange of a and lo meson, respectively. Although in chiral Lagrangians one cannot attach mean fields to Goldstone 
bosons, in the sense of Ref. , one might try to describe the effects of explicit chiral symmetry breaking by a scalar 
field a. The Lagrangian is. 


Ck = d^Kd^K - {m\ - g^KmK(j)KK + dt K, ( 22 ) 

where a and wq are the scalar and the (time-component of) vector fields, respectively, g^x and g^x are the coupling 
constants between the kaon and the scalar and the vector fields, respectively. The energies of kaon and antikaon in 
nuclear medium are then given by 


Lox = [m\ -I- - gaKmxcr + {guKtoo)'^] + gujKUJo, 


(23) 


ujx = [rrix + - gaKmxcr + {gujK^^o f] - gu^KOJo- 


(24) 


The results in the meson-exchange picture also have model dependence. First, they depend on the strength of the 
scalar and vector fields, which might be constrained from the nuclear matter properties. Second, they depend on the 
coupling constants g^x and guix- In the simple quark model, g^K/guN ~ 1/3, since there is only one light quark 
in kaon. The relation for scalar couplings is more subtle, since there is no fundamental scalar meson; it represent 
correlated s-wave two-pion exchange. 

By analysing kaon-nucleon scattering in the meson-exchange model, and by fitting to experimental data, one may 
be able to determine these coupling constants. This is the approach adopted in Refs. [ ^2|J^ by the Jiilich group. The 
coupling constants so determined, however, cannot be used directly in Eqs. ( p^ and (|24|), as these expressions are 


obtained in the simple mean-field approximation, while in Refs. 331, more much complete diagram sets and infinite 


summation of ladder diagrams are involved. This is analogous to the case of nuclear matter saturation, namely, one 
cannot apply directly the coupling constants in the Bonn model to the Walecka-type mean-field model for nuclear 
matter saturation. 

Starting from the Bonn model for nucleon-nucleon interaction, one can carry out Dirac-Brueckner-Hartree-Fock 
calculation to study nuclear matter properties. In the same sense, one can study kaon properties in nuclear matter 
starting from the Jiilich kaon-nucleon (and antikaon-nucleon) interactions, by solving in-medium scattering equation 
(the so-called G-matrix). This kind of approach should provide a consistent way of study hadron systems containing 
strangeness degrees of freedom. 


III. KAON PRODUCTION IN HEAVY-ION COLLISIONS 

One of the most important ingredients in the transport model study of particle production in heavy-ion collisions 
is the elementary particle production cross sections in hadron-hadron interactions. At SIS energies, the colliding 
system consists mainly of nucleons, delta resonances, and pions. We need thus kaon and antikaon production cross 
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sections from nucleon-nucleon {NN), nucleon-delta (NA), delta-delta (AA), pion-nucleon {ttN), and pion-delta (ttA) 
interactions. In addition, the antikaon can also be produced from strangeness-exchange processes such as ttY —>■ KN. 
In the next two subsections, we will discuss these cross sections. Because of the lack of experimental data, especially 
near production threshold that are important for heavy-ion collisions at SIS energies, we have to adopt the strategy that 
combines the parameterization of experimental data with some theoretical investigations and reasonable assumptions 
and prescriptions. 


A. Elementary kaon production cross sections 

In this subsection we discuss kaon production cross sections in pion-baryon and baryon-baryon collisions. 


1. kaon production in pion-baryon interactions 


Two major processes for kaon production in ttN interaction, namely, ttN —>■ AK and ttN —>■ YK, have actually been 
measured quite extensively in the literature |23|] . Cugnon and Lombard proposed the following parameterizations for 
experimental data in some specific channels |51|, 


^ 7T~p- 


— 


9.89(v^ — ^sq) mb if = toa -f rriK < ^/s < 1.7 GeV 
mb if Vi > 1.7 GeV 

Vs-i.e V 


(25) 




^E+K+ 


f 3.21(Vi — Vso) if Vso = WE -I- rriK < \/s < 1.9 GeV 
~ 1 7 wb if \/i > 1.9 GeV 


if Vs > 1.9 GeV 


(26) 


K+ + = 0-25 [l.O — 0.75(Vi ~ 1.682)] mb y/s > 1.682 GeV. 


(27) 


Furthermore, under some assumptions, they obtained the isospin averaged cross sections 

1 

C^nN^AK = -^(^TT-p^AKO, 


(28) 


CTt^N^SK — 2 [^r-p- 


-K+ + O'tt-p^EOA’O + ■ 


(29) 


These parameterizations have been used by several groups in the study of kaon production in heavy-ion collisions at 
1-2 AGeV region g. 

Recently, Tsushima et al. studied kaon production in ttN and ttA interactions in the resonance model. The 
kaon is produced mainly in the s-channel processes from the decay of baryon resonances such as V(1650), V(1710), 
V(1720), and A(1920). The A(1920) in their model is treated as an effective description of the contributions from 
six A-resonances between 1.9 and 1.94 GeV that have decay possibilities into YK final state. Overall, very good 
agreement with experimental data has been achieved. One of the advantages of a model calculation is that the 
isospin average can be done explicitly. In addition to kaon production in the ttN interaction, these authors have also 
calculated kaon production cross sections in ttA interactions, within the same model. Generally, these cross sections 
are substantially smaller than those in the ttN interaction, because the branching ratios of these resonances into ttA 
are smaller than those into ttN. In Fig. we show the isospin-averaged cross sections from Ref. [Q, which will be 
used in this study. For the ttN interaction we also show those from Gugnon and Lombard . 

In addition to processes without the pion in the final state, the kaon can also be produced together with one or more 
pions. At SIS energies, these processes might not be very important, since they have a higher production threshold 
and are thus suppressed as compared to those without pions in the final states. They, nevertheless, may not be totally 
negligible, and become increasingly important as beam energy increases. Here we will consider processes with one 
pion in the final state. 

For the AttK channel, there are three sets of data available, and they can be approximately fitted by the following 
expression. 


^7r~ p — ^Att^K^ 


' p—^A-K~ KA 


' ^'7r+p^A-K+ K+ 




(30) 
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where -^5 is the available energy and = m\ + mTr+mK- The comparison of this parameterization with experimen¬ 
tal data is shown in Fig. IJ. Assuming that all the other channels have the similar cross sections, the isospin-averaged 
cross section is then obtained, 




40 




4.24 


(31) 


Similarly, the six sets of available experimental data for J^ttK final state can be approximately fitted by the same 
expression 




19 




(32) 


where ,/so = ms + m-^ + mx- The comparison of this parameterization with data is shown in Fig. ^ Under the 
same assumption that all other channels have the similar cross sections, we get the isospin-averaged cross section 





(33) 


2. kaon production in baryon-baryon collisions 

The situation of experimental measurement of kaon production in nucleon-nucleon collisions is much worse than 
for pion-nucleon interactions. Until very recently, there are basically no data available near the threshold. The first 
detailed analysis of the experimental data was carried out by Randrup and Ko They proposed the following 
parameterizations for the available data 


^pp^pT,^ K+ 


^pp—^pAK+ 




= 0.024^111^ mb (34) 

rriK 

= 0.012^111^ mb, (35) 


rriK 


where Pmax is the maximum momentum of kaon given by 


Pmax — 



(ruAT -I- my + mKY)) (s — {my + my — mK^), 


(36) 


with the total energy and my the mass of hyperon in the final state. 

Based on experimental data and isospin analysis, Randrup and Ko also obtained the following isospin-averaged 
cross sections 


aNN^NYK = 0.072^ mb, (37) 

rriK 

where Y is either lambda or sigma hyperon. Furthermore, they have analysed the experimental data for processes 
with one pion in the final state, and by using the detailed-balance relation, they obtained the following expressions 
for kaon production in A^A and A A collisions, 

(JNA^NYK = 0.054^^ mb, (38) 

rriK 

ctaa^nyk = 0.036^^ mb. (39) 

mK 

Thus at the same center-of-mass energy, kaon production cross section in the A A collision is 3/4 of that in the NN 
collision, while that in the AA collision is half of that in the NN collision. 

Later on, Schiirmann and Zwermann | |^ proposed another parameterization of experimental data which treated 
more accurately the threshold behavior by using a quartic Pmax dependence rather the linear dependence of Randrup 
and Ko, 


(Jpp^K+X = “b. 


(40) 
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This parameterization was fitted to the inclusive K'^ production cross section in proton-proton (jjp) collisions. But it 
has often been identified with the lowest threshold process NN ^ NAK^. 

Recently, some of the much needed experimental data for kaon production in pp collisions near threshold have 
become available from the COSY-11 collaboration They have measured the cross section (Jpp_p\K+ at 2 MeV 

above the threshold, and obtained a cross section of about 8 nb. Based on these as well as early experimental data, 
Gassing et al. proposed new parameterizations for exclusive kaon production cross sections in pp collisions. 


^pp^pAK+ 

= 0.732 ( 

5 > 


(41) 

K+ 

= 0.275 ( 

i-f^' 

S > 


(42) 

^pp—^pH+K^ 

= 0.338 ( 

5 > 


(43) 


where are corresponding threshold. The isospin factors and the scaling factor for TV A and A A in these new 
parameterization are the same as those in Randrup-Ko parameterization. 

There have been also several theoretical calculations of kaon production in baryon-baryon interactions, based mostly 
on one-boson-exchange model p9|-|6^ . We discuss here mainly the model of Ref. , which includes one-pion and 

one-kaon exchanges. By adjusting two cut-off parameters Ajr and very good agreement with experimental data 
on the kaon production cross section in NN interactions has been achieved. In Fig. @, the results for pp pAK'^ are 
compared with experimental data and various parameterizations. The open circles are early data from the compilation 
of Baldini et al, while the solid circle is the experimental data available only very recently from the COSY-11 
collaboration that measured kaon production in proton-proton interactions at 2 MeV above the threshold . It is 
seen that our model provides a good description for both the old and new data. The parameterization of Randrup and 
Ko is shown in the figure by the dashed curve. It describes the older experimental data (open circles) reasonably 
well. It however significantly overestimates the newest data point at 2 MeV above the threshold (solid circle). The 
assumed linear dependence on Pmax (the maximum momentum of kaon) is responsible for this incorrect threshold 
behaviour. The dotted curve in the figure gives the recent parameterization of Gassing et al. |^], who has specifically 
tried to fit the newest data from the COSY-11 collaboration. 

The extension of this model to kaon production in nucleon-delta (TVA) and delta-delta (AA) collisions is straight¬ 
forward, except for one complication arising from the fact that in the AttTV vertex energy momentum conservation 
allows the pion to go on-shell. A pole develops in the pion propagator at part of the kinematically allowed region of 
the phase space. As a result the cross section becomes singular. This kind of singularity also appears in TVA ^ NNt] 
[Q, NA NN(j) 1^, and irp irp via pion exchange. The same kind of singularity occurs in the process 
—> ePW''^, in which the p can decay to and the exchanged neutrino goes on-shell in the diagram, leading 
to a singular cross section . 

In Ref. 1^, we applied the Peierls method |^], which takes care of the finite lifetime of the incoming A resonance, 
to regulate the singularity. The singularity is removed by the delta width gained through the energy-momentum 
conservation relation. This width is density independent and the resulting cross section does not diverge as the 
density goes to zero. The Peierls method has also been used in Ref. |^] to remove the singularity associated with 
muon collision. 

The results for isospin-averaged kaon production cross sections are shown in Fig. ^ together with those from 
the parameterization of Randrup and Ko | |^ , with the upper curves represent the ANK final state, and the lower 
curves the UNK final state. For the NN interaction, the Randrup-Ko parameterization is larger than our results, 
since in the former only one-pion-exchange was considered which leads to a large isospin-average factor. For TV A and 
AA interactions, our results are generally larger than the Randrup-Ko parameterization, except near the threshold. 
To apply these cross sections in the transport model, we have parameterized our theoretical results in terms of the 
following expression 


<^BB^NYK 


a(ys-v^)^ , 

■;-r-F- --- mb 

5-1- (Vs - 


(44) 


The fitted parameters a, b and x for six channels are listed in Table 1. These parameterizations are shown in Fig. 
by dashed curves, which are seen to reproduce the solid curves quite accurately. 


Table 1 Fitted parameters in Eq. (^^. 


9 



NN NAK NN NY,K NA NAK NA NYsK AA NAK AA NYK 


a 

0.0865 

0.1499 

0.1397 

0.3221 

0.0361 

0.0965 

b 

0.0345 

0.167 

0.0152 

0.107 

0.0137 

0.014 

X 

2.0 

2.4 

2.3 

2.3 

2.9 

2.3 


In addition to processes with no pions in final states, kaons can also be produced together with one or more pions 
in the final states. We have analysed the available experimental data, and propose the following parameterizations 
from processes with one and two pions in final states. For pp —f NAttK, the three sets of experimental data can be 
fitted by the same expression 


^pp^Ap 7 T+ ~ ^pp—*Ap'iT'^K+ ~ >An 7 r+iC+ 


0.084(Vi - 

0.45 + (v/i-“ ■ 


(45) 


We assume that other proton-proton (pp) and neutron-neutron (nn) channels have the same cross section. The only 
set of data that is available for neutron-proton (np) channels indicates a cross section which is about half of the pp 
one, and we assume all other np channels have the same cross section as this np channel. Under these assumption we 
get the isospin-averaged cross section 


C^NN^NA-nK 


0.45 -I- (a/s - a/so)^ 


The comparison of these parameterizations with experimental data is shown in Fig. 

For pp NYttK, on the average, the five sets of available data can be fitted by the following expression 


Ir , ,1^, 0.04(Vi-v^)2 ^ 

5 [^pp^p^+nOKO + <ypp^p^+^-K+ + •••]- + (v^ - 75^)2 


(46) 


(47) 


We assume that all the other pp and nn channels have the similar cross sections. For np channels the three sets of 
available experimental data can be fitted, on the average, by the following expression, 

■032(V^-v^)^ 

35 + (Vi 






If A 0 

^ [^np— ^pl2^'K~K^ “b ^np —*pS‘^7r“/C+ ~b ^np —^pcr“7r+iC^J ~ ^ ■ 


We assume that all the other np channels have the similar cross sections. We thus get the isospin-averaged cross 
section. 


<yNN^NY.-nK 


o.i6(Vi-yi^)^ o.i6(Vi-yi^)^ 
0.75 + (yi - y5F)2 1.35 + (yi - v^)2 


(49) 


The comparisons of these parameterizations with experimental data are shown in Fig. 

Finally, we discuss cross sections with two pions in the final state. For pp NAttttK, the available three sets of 
experimental data can be fitted, on the average, by the following expression, 


1 

3 


[^pp—*pA7r+7r‘^iC*^ “b ^pp—^ pAtt+tt iC+ “b O'pp —,pA'7r+7r+iC‘’] 


o.o8(yi - v ^)2 

2.25 + (yi-y5^)2 “ • 


(50) 


We assume that other pp and nn channels have the same cross section. The only set of data that is available for np 
channels indicates a cross section which is about half of the pp one, and we assume all other np channels have the 
same cross section as this np channels. Under these assumptions we get the isospin-averaged cross section 

o.32(yi - yi^)2 ^ 

crNN^NA^.K ^ 2.25 + (yi-y5F)2 

The comparison of these parameterizations with experimental data is shown in Fig. |^. 

A similar procedure is carried out for NN —> NYttttK. The six sets of data for pp channels can be fitted, on the 
average, by the following expression. 


ir_ ,1.. o.o5(yi-yi^)2 

g ypp^pE+TT+TT K° + Cpp^pS+TrOTT K+ + -I 2 25 + (yi — Y^)^ 


(52) 
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Again, the np channels have a cross section which is about half of that for pp channels. Under the same assumptions 
as before, the isospin-averaged cross section is found to be 




0.55(Vs - 

2.25 + (yi-ys^)2 


(53) 


The comparisons of these parameterizations with experimental data are shown in Fig. 0 

Combining our theoretical results for pp —>■ NYK'^ with our parameterizations for channels with one and two 
pions in final state, we can calculate the inclusive K'^ production cross section in pp collision which is compared with 
experimental data in Fig. 12, In the figure, the dotted, short-dashed, and long-dashed curves are for final states with 
zero, one and two pions, respectively, while the solid line is the sum of the three contributions. It is seen that our 
results are in good agreement with the data. In the figure we also show the Schiirmann-Zwermann parameterization 
which was originally proposed for the inclusive production cross section in pp interaction. 


B. Elementary antikaon production cross section 

Similar to the kaon, at SIS energies, the antikaon can be produced from ttTV, ttA, NN^ NA, and AA interactions. 
It can also be produced from the strangeness-exchange processes such as ttY KN. Although the abundance of 
hyperons in heavy-ion collisions at SIS energies is small, these strangeness-exchange processes are quite important, 
because of a large cross section. 


1. antikaon production in pion-baryon collisions 


The cross sections for ttN NKK have been analysed recently by Sibirtsev et al. using a boson-exchange 
model. The authors found that the ratios between experimental cross sections of different charge channels can be well 
accounted for by this model. The isospin-averaged cross section from their model is 


^ttN^NKK — SCTti— if-, 

with the latter taken from the parameterization of the experimental data. 


(54) 


a^-p^pKOK- = 1-121 (^1 - mb. 


(55) 


where = rriN + 2m k- The comparisons of this parameterization, after including correct isospin factor for different 
charge channels, with the experimental data are shown in Fig. |^. 

There are also processes with one pion in final states. The available experimental data indicate charge independence 
within a factor of two, and can be parameterized by 


^7r — p — ^p-TT® if 0 if - ~ — p — *p7T ~ K+ K~ 


85.0 



^£o\ 

V s ) 

^ s ) 


mb, 


(56) 


where y/so = mN + rrin + 2mK- The comparisons of this parameterization with experimental data are shown in Fig. 
[l^ . The isospin-averaged cross section is 

cr^N^N^KK = 510.0 (l - (y) mb. (57) 

There has been no theoretical study on antikaon production in pion-delta collisions. We assume that these cross 
sections, after isospin averaging, are the same as those for ttN interactions at the same center-of-mass energy. 


8. antikaon production cross sections in baryon-baryon collisions 

Again, the situation of experimental data for antikaon production in NN collisions is worse than that for ttN 
interactions. Zwermann and Schiirmann [T^ analysed the available experimental data and, under the assumption 
that the cross section is charge independent, proposed the following parameterization for the isospin-averaged cross 
section. 
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^NN^NNKK — O.OSpmax mb, 


(58) 


where Pmax is the maximum momentum of antikaon given by 


Pmax — 


2^5 


\J{s- 4(TOAr + mKY){s - Am%). 


(59) 


Recently, Sibirtsev et al. analysed this process in a boson-exchange model, using as inputs the antikaon pro¬ 
duction cross sections from pion-nucleon interactions outlined above. Their results fit the experimental data for 
pp ppK^K^ quite well, but seem to underestimate the data for pp npK^K^. We use in this work the following 
parameterization for this process 


/ 5n \ ^ / Sn\ 0-31 

<^pp^npK+K0 = 0-19 - — j ( —j mb. (60) 

Experimental data indicate that 

^pp^ppK^K^ ~ ^np^ppK^K~ ~ ~^^pp^npK+■ (^1) 

The comparision of these parameterization with experimental data are shown in Fig. M. Note that at 2 MeV above 
the threshold, our parameterization gives a cross section of about 0.09 nb for pp ppI^K^. This is very close to the 
preliminary data of about 0.1 nb from the COSY-11 collaboration for pp ppK^K~ Both the and K~ 
data near the thresholds from the COSY-11 collaboration are indeed very helpful in constraining the inputs of and 
reducing the uncertainties in the transport models. With some further assumptions about the relationship between 
different charge channels, we obtain the following isospin-averaged cross section 

5 

^NN^NNKK = pp^npK+K° ■ ( 62 ) 

The comparision of this with the parameterization of Zwermann and Schiirmann is shown in Fig. |^. Apparently, 
the latter is much larger than ours near the production threshold. Again the linear Pmax dependence assumed in the 
Zwermann-Schiirmann parameterization is responsible for the incorrect threshold behavior. It should be mentioned 
that the Zwermann-Schiirmann parameterization is for the isospin-averaged cross section Unn^nnkR- In 
it was incorrectly compared with the inclusive K~ production cross section in pp collision. 

For processes with one pion in the final state, experimental data are quite scarce. We propose the following 
parameterization for pp —> pmr'^K^K^, 

(^pp^pn^+KoRo « 0.945 (^1 - mb. (63) 

The available data for other channels indicate that 


^pp^ppn'^ K^K~ ‘ 


^ 1 

■' '^^pp^pmT+■! 

(64) 


^np^ppTT- ^ 

^ 1 

'^^pp—*pn7T+K^K^’ 

(65) 


The comparisons of these parameterizations with experimental data are shown in Fig. HI- With some further 
assumption about the relationship between different charge channels, we obtain the following isospin-averaged cross 
section 


^NN^NNttRR 



( 66 ) 


Finally, on the average, the available experimental data for five channels with two pions in the final state can be 
fitted by the following expression, 


1 

5 


— *pp7T'^ 7T^ R^ R “f ^pp —^priTT+TT R'^ R^ -f * * ’] ~ 0.316 



(67) 


The comparison of this parameterization with experimental data is shown in Fig. E- The isospin-averaged cross 
section, assuming that all the other channels have the same cross section, is given by 
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^NN^NN-irnKK 




mb. 


(68) 


Summing up the cross sections with zero, one and two pions in the final state, we obtain the inclusive K~ production 
cross section in pp collision, which is compared with experimental data in Fig. |^. The dotted, short-dashed, and 
long-dashed curves correspond to channels with zero, one and two pions, respectively. The sum of these contribution, 
as shown in the figure by solid curve, agrees rather well with the data. 

It is also useful to compare inclusive K'^ and K~ production cross sections in pp collisions. This is done in Fig. 

It is seen that at the same \/s — y/so {y/so = mjsf + toa -I- ttik for case, and = 2mN + 2mK for K~ case), 
the production cross section is one to two orders of magnitude larger than that of K~ in pp collisions. 


3. antikaon production in pion-hyperon collisions 


As mentioned, in heavy-ion collisions, the antikaon can also be produced from strange-exchange processes such as 
ttA —> KN and ttE ^ KN. These cross sections are not known empirically, but they can be obtained from the inverse 
processes, which have been studied in great detail experimentally, by using the detailed balance relation. At low 
energies, the experimental data on KN —> ttY can be analysed in terms of the K-matrix for three coupled channels 
KN, Att, and Ett |49 Based on this model, Ko obtained the cross sections for ttA —s- KN and ttE ^ KN which 
can be as large as a few mb [Q. 

In Ref. Cugnon et al. proposed the following parameterizations for the experimental data. 


J 50p2 - 67p + 24 mb if 0.2 < p < 0.9 GeV 
ctk-p^At^o - j 3 p- 2.6 jf 0.9 < p < 10 GeV/c ’ 

(69) 

= 0 . 6 p ^■®, 

(70) 

J 1.2p“^-^ mb if 0.5 < p < 1 GeV/c 

CTi 7 -n^E 0 ^- - 1 1 2p-2-3 mb if 1 < p < 6 GeV/c ’ 

(71) 


where p is antikaon momentum in the laboratory frame. These parameterization do not describe the resonance 
structure in the data very well, as will be shown in Fig. We will use the following parameterizations instead, 

p < 0.6 GeV 

0.6 < p < 1 GeV/c , (72) 

p > I GeV/c 


if p < 0.345 GeV 

« cTic-n^EOTT- = { 0.0138/((p- 0.385)2 q.OOI?) mb if 0.345 < p < 0.425 GeV/c . (73) 

[ 0.7p“2 09 mb if p> 0.425 GeV/c 

The comparisons of these parameterizations with experimental data are shown in Fig. The dotted lines give the 
parameterizations of Cugnon et al. |^. According to Ref. |Q, the isospin-averaged cross sections for the inverse 
processes are then 


O’A'-P^AttO 


I.205p-i-428 mb 
= 3.5p° ®^® mb 

3.5p-3-97mb 


0.624p ^-88 mb 


^nA^KN — 2 



2 


^K~p —»^A7r° 7 


(74) 




fP^Yf , 


O'K-r 




r-), 


(75) 


where px and p^^ are, respectively, momenta of antikaon and pion in the center-of-mass frame. These cross sections 
are shown in Fig. |2^. At low invariant energies, they are in agreement with those obtained in Ref. based on the 
K-matrix method. The first peak in these cross sections comes from the threshold effect, while the second peak comes 
from the resonance structure in the KN —> ttY cross sections. 
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C. Kaon and antikaon final-state interactions 


Particles produced in elementary hadron-hadron interactions in heavy-ion collisions cannot escape the environment 
freely and be detected. Instead, they are subjected to strong final-state interactions. For the kaon, because of 
strangeness conservation, its scattering with nucleons at low energies is dominated by elastic and pion production 
processes, which do not affect its final yield but changes its momentum spectra. Kaon-nucleon scattering was been 
studied in detail experimentally |74|. Data were nicely summarized in Ref. |49|. 

The experimental data for elastic kaon-nucleon scattering can be parameterized by 


<^KN^KN = 


10.5-f 19.7(Vs-mb 


if < 0.2 GeV 


5.76(Vs - V^)/(0.056-k (v^- V^)2) mb if ys-v^>0.2GeV 


(76) 


where = uin + rriK- 

Similarly, the one-pion {KN KNtt) and two-pion production {KN KNtttt) cross sections can be parameter¬ 
ized as 


OKN^KN. = 1.557(1 - f£)3.28(f0 )4.35 ( 77 ) 

s s 

with = niN -\- rriK + and 

(^KN^KNtttt = 12.87(1-2.^6.38^_0^4.9 

s s 


with = rriN + ttik + 2m,r. Tbe comparison of these parameterizations with experimental data is shown in Fig. 

The final-state interaction for the antikaon is much stronger. As mentioned, antikaons can be destroyed in the 
strangeness-exchange processes. They also undergo elastic scattering. We parameterize the experimental data for 
K~p total, elastic, and charge-exchange cross sections. 


^K~p—>total — 


23.5p-i-°4 mb if 

0.504/((p-0.39)2-p 0.0056) mb 
181.9(p-0.75)2-p 34.0 mb 
55.2p-i'85 mb if 


p < 0.35 GeV 
0.35 <p< 0.46 GeV/c 
0.46 <p< 1.05 GeV/c ’ 
p > 1.05 GeV/c 


(79) 


J 11.2p-°'986 mb if p < 0.7 GeV 

I 5.0/((p-0.95)2-f 0.25) mb p > 0.7 GeV/c ’ 


(80) 


1.813p mb if p < 0.35 GeV 

0.0192/((p-0.39)2 -I- 0.0016) mb 0.35 < p < 0.43 GeV/c . (81) 

15.9/((p - 0.9)2 2.65) mb if p > 0.43 GeV/c 


where p is the momentum of the antikaon in the laboratory frame. The comparisions of these parameterizations with 
the experimental data are shown in Fig. |^. At low momenta, our parameterizations for the total and elastic cross 
sections are very similar to those proposed in Ref. ||4^ ]. For the K~n interaction, experimental data are available only 
at relatively large momenta. These data indicate the following approximate relations, 


^K~ n^total 


(82) 

^K-n^K-n " 

0.9(7/^ —p— >K~p’ 

(83) 


The comparision of these relations with the experimental data is shown in Fig. ^ 

We treat the charge-exchange process K~p <-> K^n as ‘elastic’, since we do not have explicit isospin degrees of 
freedom in our transport model. Under these approximations, the isospin-averaged KN elastic cross section is given 
by 


^KN^KN ~ 0.5(1.9crx-p^iC-p + ^K--p^K0n)y (84) 

and the antikaon absorption cross section is given by 

^KN^absorption ~ ^K~p^total ^KN^KN' (85) 

Both the elastic and absorption cross sections increase rapidly with decreasing antikaon momenta. This will have 
strong effects on the final K~ momentum spectra in heavy-ion collisions. 
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D. Some Discussions 


1. Inclusive versus exclusive cross sections: As mentioned, kaon (and antikaon) can be produced in final states with 
zero, one, and more pions. To calculate their total production probability, inclusive cross sections such as pp K'^X 
and pp K~X should be used. On the other hand, to treat their final-state interactions and to calculate their 
spectra, exclusive cross sections with specific final states, such as pp pAK'^ and pp pAirK'^ are needed. For NN 
collisions, we have shown that by including exclusive processes with zero, one, and two pions in the final states, the 
inclusive cross sections for both AT'*' and K~ production with invariant energy up to about 5 GeV can be saturated. 
For ttN interaction, we have included exclusive processes with zero and one pion in final states. They are expected to 
saturate the inclusive cross sections with invariant energy up to about 3 GeV. Thus, as far as the heavy-ion collisions 
at SIS energies are concerned, we have included sufficient exclusive processes for both the NN and irN interactions. 

2. Total versus differential cross sections: So far we have discussed only total cross sections for kaon and antikaon 
production in a specific channel, which give their production probabilities. To calculate their spectra and to treat their 
final-state interactions, we need also to determine their momenta. We need thus also differential cross sections, such 
as momentum spectra and angular distributions. There are far less experimental data on differential cross sections 
than on total cross sections. Model calculations trying to describe both the differential and total cross sections are 
also far more involved than the one-boson-exchange model of Ref. |61| that attempted at total cross sections only. 

For two-body final states, such as irN YK, the magnitude of the kaon momentum in the ttN center-of-mass frame 
is fixed. We thus need only angular distribution to determine the direction of the momentum. In all the transport 
model calculations, this is assumed to be isotropic in the ttN center-of-mass frame. We will also use this assumption in 
this work. It is expected that, because of strong final-state interaction, the use of an anisotropic angular distribution 
will not have significant effects on the final kaon momentum spectra. To see more quantitatively the effects of angular 
distribution, we will take ttN —> YK as an example, where limited experimental data are available [fz^ . In Fig. 
we show the angular distribution for Tr~p AK^ at = 1.742 GeV. The open circles are experimental data, and 
the solid curve is our parameterization. 


^ = O.OI -I- O.OI75(cos0 -I- 0.93)^ mb/sr. 
as 2 


( 86 ) 


We neglect the beam energy dependence of the angular distribution. The results for kaon final spectra using this 
angular distribution will be compared with those using isotropic assumption. 

For three-body final states, such as pp —>■ pAK'^, only the maximum momentum of the kaon is fixed by energy- 
momentum conservation. Usually one first determine the magnitude of the kaon momentum in the pp center-of-mass 
frame by using some sort of momentum spectra. Randrup and Ko proposed the following momentum spectra, which 
has been frequently used in transport models. 


da 

dp 



(87) 


where Pmax is the maximum momentum of the kaon in the NN center-of-mass frame. This parameterization is seen 
to describe the experimental data quite well, as can be seen from Fig. which will also be used in the present 
work. To see the effects of elementary kaon momentum spectra on the final kaon momentum distribution in heavy-ion 
collisions, we will also use a somewhat different parameterization, which also provides a reasonable description of the 
data. 


da 

dp 



( 88 ) 


This parameterization is also compared with the data in Fig. 

3. In-medium versus free-space cross sections: As cab seen from last subsections, we have parameterized all the 
elementary production cross sections in terms of ^/s and where y/s is the available energy and is the 

threshold. In nuclear medium, kaon and antikaon masses are modified, so are their production thresholds. We will 
then use \/s^, which are calculated with effective masses, in evaluating the in-medium production cross sections. 
This amounts to the change of threshold, or approximately, to the change of available phase space. This kind of 
treatment of medium effects is certainly incomplete. In nuclear medium, not only hadron masses, but also coupling 
constants, cut-off parameters, and even the structure of the cross sections might change. Nevertheless, the change of 
the thresholds, or phase space, is the most apparent and the simplest to implement. To study other medium effects, 
one needs a complete theory for the elementary cross sections. 
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Since the K~^ mass increases and the K~ mass decrease in nuclear medium, the K'^ production cross sections are 
suppressed while those of K~ are enhanced in heavy-ion collisions, when medium effects on their masses are included. 
This leads to an increase of K~ /ratio with nuclear density, as illustrated schematically in Ref. [ 0 - This increase 
is partly responsible for the observed K~/K~^ ratio by the Kaos collaboration ||2^. We will come back to this point 
later. 

In addition to the change in the production cross sections, the medium effects on kaon and antikaon also affect their 
momentum spectra, when they propagate in the mean-field potentials. The Hamilton equations of motion for kaon 
and antikaon are very similar to those for nucleons. 


dr 

dt 


k dk 

T bkPN dt 


— '^xUk^Kj 


(89) 


where minus sign corresponds to kaon, and plus sign to antikaon. It is clearly that the K'^ momentum increases and 
that of K~ decreases when they propagate in their respective mean field potentials. This is affect significantly their 
momentum spectra, and especially the momentum spectra of their ratio. 


E. Relativistic transport model 


Heavy-ion collisions involve very complicated nonequilibrium dynamics. In order to extract useful information on, 
e.g., nuclear equation of state and hadronic properties in dense matter, it is necessary to use transport models. At 
SIS energies, both mean-field and two-body scattering play important roles in the dynamical evolution of the system 
and need to be included. Over more than ten years, Boltzmann-Uehling-Uhlenbeck (BUU) and similar models based 
on either Skyrme-type or Walecka-type effective nucleon-nucleon interactions have been developed |19 In 

this work we used the relativistic transport model (RVUU) similar to that developed in Ref. and used in Refs. 
[^,^. Instead of the usual linear and non-linear tr-w as in Refs. we base our model on an effective chiral 

Lagrangian recently developed by Furnstahl, Tang, and Serot [^, which will also be used in our studies of neutron 
star properties. 

The effective chiral Lagrangian of Furnstahl, Tang and Serot is derived using the dimensional analysis, naturalness 
arguments, and provide a very good description of nuclear matter and finite nuclei. In the mean-field approximation, 
the energy density for the general case of asymmetric nuclear matter is given by 




e = 


( 2^)3 


dk\/k‘^ 






{2^f 


dk\lk? + m*^ 


-\- W p ~\-R Pn) 2C^ 






4C| 


1 - 


S' 


4/rf r 


5‘” (‘ -1 


_ - ^ — 

24 2Cl S' 


IFA 


(90) 


In the above the first line represents the ‘kinetic energy’ of nucleons, with Kp^ and Kp^ being proton and neutron 
Fermi momenta determined from their densities, Pp and pn- The nucleon effective mass is related to its scalar 
field $ by m*j^ = mjv — 4>. In the second line, W (from uj meson exchange) and R (from p meson exchange) are the 
isospin-even and isospin-odd vector potentials, respectively. The latter vanishes for symmetric nuclear matter. The 
third line includes the self-interactions of the scalar field (proportional to the anomalous dimension d ), the vector field 
(proportional to ^), and the coupling between them (proportional to p). The parameters of this model are adjusted 
to fit the properties of nuclear matter and finite nuclei. We will use the parameter set T1 listed in table 1 of Ref. [ pl| . 
The nucleon scalar and vector potentials for symmetric nuclear matter are shown in Fig. ^ as a function of density. 
Note that because of self-interaction of the vector field and the coupling between the scalar and vector fields in this 
model, the vector potential tends to saturate at high densities. In the usual linear Walecka model and non-linear 
(j-uj models , the vector potential increases linearly with density. In Ref. jsj] we have shown that the inclusion of 
higher-order terms for vector field is instructive in reproducing correctly both the nucleon flow and dilepton data. 

From the energy density, we can derive a relativistic transport model for heavy-ion collisions. At SIS energies, the 
colliding system consists mainly of nucleons, delta resonances, and pions. While medium effects on pions are neglected 
as in most transport models, nucleons and delta resonances propagate in a common mean-field potential according to 
Hamilton equation of motion, 


dx p* dp 

dt E* ’ dt 




(91) 
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where E* = + In our transport model we have neglected the explicit isospin degrees of freedom and 

thus the contributions from p meson exchange. In addition to propagations in their mean field potentials, we include 
typical two-body scattering processes such as BB ^ BB, NN <-> NA and A Ntt. Cugnon parameterizations and 
proper detailed-balance prescriptions are used for describing these reactions . 


F. Proton and pion observables 

To show that our transport model describes reasonably the nucleon, delta and pion dynamics, we compare our 
predictions for proton and pion observables with recent data from the FOPI collaboration. In Fig. |2^ we show 
proton transverse mass spectra in the mid-rapidity region in central Ni-|-Ni collisions at 1.06, 1.45, and 1.93 AGeV. 
The impact parameter 6 < 2 fm is chosen so that fair comparison with FOPI data [^, which has a geometry cross 
section of about 100 mb, can be made. It is seen that our results for 1.06 AGeV agree very well with FOPI data 
[^ . Our spectra can be fitted by exponential function exp(—wt/T). The slope parameter T is found to be about 95, 
115, and 140 MeV at 1.06, 1.45 and 1.93 AGeV, respectively. These are also in good agreement with what have been 
extracted from FOPI data 

In Fig. ^ we show our results for 7r“ transverse mass spectra in the same reactions as in Fig. ^ Again, our 
results are in very good agreement with FOPI data for central Ni-|-Ni collisions at 1.06 AGeV. The slope parameters 
of our pion spectra, fitted from mt — m=0.1 to 0.5 GeV, are about 70, 85, and, 95 MeV for 1.06, 1.45, and 1.93 AGeV, 
respectively. The difference between proton and pion slope parameters can be explained by the collective transverse 
flow which has stronger effects on the heavier protons than pions . 

In Fig. ^ we compare our predictions for the proton rapidity distribution in central Ni-|-Ni collisions with FOPI 
data. Here the solid circles are experimentally measured data, while the open circles are obtained from the measured 
data by reflecting with respect to the mid-rapidity. As in the FOPI data, in this paper y is defined as the rapidity of a 
particle in the nucleus-nucleus center-of-mass frame, normalized by the beam rapidity. Overall very good agreement 
with the data have been achieved. Similar comparisons for 7r“ rapidity distributions are shown in Fig. Again, 
the agreement with the data is fairly good. 


G. Kaon and antikaon production in heavy-ion collisions 


Since the proposal of Aichelin and Ko that subthreshold yield may be sensitive to nuclear equation of state at 
high densities [ p^ , many works have been done concerning the subthreshold kaon production in heavy-ion collisions 
using variou s transpo rt models [p^p0|j^j8^ -p^ . There have been also several investigations on subthreshold antikaon 
production 


In this paper, we will concentrate on comparisons with recent data from the KaoS collab- 
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oration for Ni-|-Ni collisions |22| using the new sets of elementary cross sections as outlined in subsections A-C. We 
consider two scenarios, namely, with and without kaon medium effects. As mentioned, we use ax ~ 0.22 GeV^fm^ 
for K'^. For K~, we adjust such that we achieve a good fit to the experimental K~ spectra. We find ~ 0.45 
GeV^fm^. With this value and the nuclear scalar density from the model of Furnstahl, Tang and Serot (the magni¬ 
tude of the scalar attractive depends not only on the value of k but also on pg, which in turn is model dependent 
1^), we find that the K~ feels an attractive potential of about 110 MeV at normal nuclear matter density, which is 
somewhat smaller than that extracted from kaonic atoms | p^ . The K~ production is not, however, sensitive to the 
lower densities as those probed in kaonic atoms. 


1. Sensitivity to differential cross sections 

Before presenting our main results and comparing with experimental data, we discuss in this subsection the sen¬ 
sitivity of final spectra to the differential cross sections. We calculate transverse mass spectra and rapidity 
distribution in Ni-|-Ni collisions at 1.8 AGeV and b=0 fm. For the sensitivity to the angular distribution we take 
ttN —^YK as an example. In the first case, we assume that kaons are emitted isotropically in the ttN center-of-mass 
frame, while in the second case, we use the angular distribution as given by Eq. The results are shown in Fig. 

As can be seen, because of rescattering, the final kaon spectra is not very sensitive to the elementary angular 
distribution. For the sensitivity to the momentum spectra, we take W —> NYK as an example. In the first case 
we use the Randrup-Ko parameterization, Eq. and in the second case we use Eq. The results are shown in 
Eig. Again, as far as the elementary momentum spectra reproduce the experimental data for pp reasonably, the 
resulting K~^ spectra in heavy-ion collisions are very similar. 
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2. and K kinetic energy spectra 


The results for kinetic energy spectra in Ni+Ni collisions at 0.8, 1.0 and 1.8 AGeV are shown in Figs. |^, 
and 37, respectively. The impact range is chosen to be 6 < 8 fm. The solid histogram gives the results with 


kaon medium effects, while the dotted histogram is the results without kaon medium effects. The open circles are the 
experimental data from the KaoS collaboration . It is seen that the results with kaon medium effects are in good 
agreement with the data, while those without kaon medium effects slightly overestimate the data. We note that kaon 
feels a small repulsive potential; thus the inclusion of the kaon medium effects reduces the kaon yield. The slopes of 
the kaon spectra in the two cases also differ. With a repulsive potential, kaons are accelerated during the propagation, 
leading to a larger slope parameter as compared to the case without kaon medium effects. 

The results for the K~ kinetic energy spectra are shown in Fig. ^ for Ni+Ni collision at 1.8 AGeV. It is seen that 
without medium effects, our results are about a factor 3-4 below the experimental data. With the inclusion of the 
medium effects which reduces the antikaon production threshold, the K~ yield increases by about a factor of 3 and 
our results are in good agreement with the data. This is similar to the findings of Gassing et al. Also, with the 
inclusion of kaon medium effects, the slope parameter of the K~ spectra decreases, since the propagation of antikaons 
in their attractive potential reduces their momenta. 

As mentioned earlier, the KaoS observation that the K~ yield at 1.8 AGeV is similar to K'^ yield at 1.0 AGeV is 
an indication of the kaon medium effects. This can best be seen in our results by looking at their ratio as a function 
of kinetic energy. This is done in Fig. By doing this we can also see more clearly the effects of kaon and antikaon 
mean-field potentials on the shapes of their momentum spectra. It is seen that without kaon medium effects, the 
K'^/K~ ratio decreases from about 7 at low kinetic energies to about 1 at high kinetic energies. In central collisions 
the ratio of their total yield is IK~ « 6, which is much greater than the KaoS observation of about 1. Since the 
antikaon absorption cross section by nucleons increases rapidly at low momentum, low-momentum antikaons are more 
strongly absorbed by nucleons than high-momentum ones. This makes the K'^/K~ ratio increase with decreasing 
kinetic energies. When medium effects are included, we find that the /K~ ratio is almost a constant of about 1 in 
the entire kinetic energy region, which is in good agreement with the experimental data from the KaoS collaboration 
p^ , shown in the figure by open circles. The ratio between their total yields in the central collisions is K'^/K~ k, 1.1, 
very close to the KaoS observation of about 1. As mentioned, the shapes of Ar+ and K~ spectra change in opposite 
ways in the presence of their mean-field potential. Kaons are ‘pushed’ to high momenta by the repulsive potential, 
while antikaons are ‘pulled’ to low momenta. If there were no antikaon absorption, or if the absorption cross section 
were independent of antikaon momentum, the /K~ would increase with kinetic energy because of propagation in 
the mean-field potential. A ‘flat’ IK~ spectra is thus highly nontrivial, which must be result from the combined 
effects of mean-field and energy-dependent antikaon absorption. An accurate experimental measurement of this ratio 
can thus be very useful in determining whether there are medium effects or not on kaons in nuclear matter. 


3. and K transverse mass spectra and rapidity distributions 


In Figs. ^ and |4l| we show the K~^ transverse mass spectra at mid-rapidity in central Ni+Ni collisions at 1.0 and 
1.93 AGeV, respectively. The impact parameter range is chosen to be 6 < 4 fm, comparable to the centrality selection 
of the FOPI collaboration’s measurement of kaon production in Ni+Ni collisions. At 1.0 AGeV, the slope parameter 
of the AT'*' spectrum is about 60 MeV when kaon medium effects are neglected. This increases to about 70 MeV when 
kaon medium effects are included. At 1.93 AGeV, the slope parameters are 95 and 110 MeV, respectively, without 
and with kaon medium effects. Similarly, the K~ transverse mass spectra in central Ni+Ni collisions at 1.8 AGeV are 
shown in Fig. The difference between the scenarios with and without kaon medium effects is shown to be most 
pronounced at low transverse mass. 

In Fig. El- we show the Ar+ rapidity distribution in central Ni+Ni collisions at 1.0 AGeV. The solid and dashed 
curves are obtained with and without kaon medium effects. The most significant difference between the two cases 
is seen around the mid-rapidities. The full width at half maximum (FWMH) is about 1.4 in the case without kaon 
medium effects. This increases to about 1.8 when kaon medium effects are included. The results for the K~^ rapidity 
spectra in Ni+Ni collisions at 1.93 AGeV are shown in Fig. 44. Here our results are compared with experimental 
data from the FOPI and KaoS collaborations [^. The FOPI measured data are shown in the figure by solid circles, 
while the open circles are obtained by reflecting with respect to mid-rapidity. The KaoS data point, shown by solid 
square, is obtained by multiplying the original KaoS data measured at 1.8 AGeV by a factor of 1.45 as in Ref. |^] 
to take account the beam energy dependence. The agreement with the data is better when kaon medium effects are 
included. In this case, while the total yield agrees with the data, our rapidity distribution is somewhat broader 
than what the data indicate. 
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We show in Fig. 45 the K~ rapidity distribution in central Ni+Ni collisions at 1.8 AGeV. The difference between 
the scenarios with and without kaon medium effects is most visible near mid-rapidity. In Fig. ^ we show the ratio 
of the fK~ as a function of normalized rapity y. As in Fig. ^ the K'^ is for beam energy of 1.0 AGeV, and 
the K~ is for beam energy of 1.8 AGeV. Without kaon medium effects, we find that the ratio decreases from about 
6.5 around the mid-rapidity to about 5 around the target and projectile rapidities. When kaon medium effects are 
included, the trend is just the opposite: the ratio increases from about 0.9 around mid-rapidity to about 2 at target 
and projectile rapidities. Experimental information on the rapidity dependence of the /K~ ratio should be very 
useful in distinguishing the two scenarios. 


4- excitation functions 

Finally, we show in Fig. the beam energy dependence of K'^ yield in central Ni-|-Ni collisions. What is shown 
in the figure is the ratio between K~^ yield and participant nucleon number Apart- Our results, shown by open 
squares connected by solid lines, are obtained with kaon medium effects included. The data from the FOPI and KaoS 
collaborations are shown by solid and open circles, respectively For comparison, we also plot in the figure 

the beam energy dependence of jApart ratio. From 0.8 AGeV to 1.93 AGeV, the pion yield increases by about 
a factor of two and half, indicating an almost linear dependence on Abeam- This is in agreement with experimental 
observations. The K'^ yield, on the other hand, increases by about two orders of magnitude. Its beam energy 
dependence is about The difference between the beam energy dependence of pion and kaon yield reflects the 

fact that at these energies kaon production is ‘subthreshold’ or ’near-threshold’, while pion production is well above 
the threshold. 


IV. KAON CONDENSATION IN NEUTRON STARS 


As mentioned, the medium modification of kaon properties affects not only kaon observables in heavy-ion collisions, 
it also bears important consequences in the structure and evolution of compact obj^ects in astr ophy sics, especially the 
maximum mass of neutron stars. There have been many studies on this subject P,[7|, p9|J^ , 1 PC |. In this paper, we 
emphasize the constraints on kaon in-medium properties from heavy-ion data as discussed in the last section. We shall 
concentrate on cold, neutrino-free neutron stars which are in chemical equilibrium under (3 decay processes. We will 
consider two scenarios: one excludes kaons and the other includes them. In the first case, we need to find the ground 
state for a system of nucleons, electrons and muons in chemical equilibrium, while in the second scenario, we try to 
find the ground state of a system of nucleons, electrons, muons and kaons in chemical equilibrium. The ground state 
energy density and pressure are then used in the Tolman-Oppenheimer-Volkov (TOV) equation to study neutron star 
properties 

In the absence of kaons, the chemical equilibrium conditions require that the chemical potentials should satisfy 




(92) 


and in the presence of kaons. 


ti = y-r 


lip — fie — flp — lix • 


(93) 


In both cases y is the overall charge chemical potential. The local charge neutrality can be imposed by minimizing 
the thermodynamical potential 


= En + £k + £l - y{Pp- Pk - Pe-Pn)- (94) 

For the energy density of nucleons we use the results of Furnstahl, Tang, and Serot, as given in Eq. (|^). The energy 
density of leptons is given by 


£L=£e + 'n[y-'mp)ep„ (95) 

where p is the Heaviside function {p{x) = 1 if x > 0 and p{x) = 0 if x < 0). The electron and muon energy densities 
are, respectively. 


Ee = 


47r^ ’ 


(96) 
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and 


m 


4 , 


Eu = 


Stt^ . 


i(l + 2i^)\/l + - ln(t + \/l + 


(97) 




where t = and Kp^ = fj,^ — mf^. The energy density of kaons, according to Ref. is 

£k = -/^ ysin^d + 2mpfsm‘^^ + nxp 

0 0 

- m(1 + x)psm‘^- - 2/^a^psin^-, 


(98) 


where / is the pion decay constant, x = Pp!p is proton fraction, and the chiral angle 0 is defined in terms of kaon 
amplitude vk. 


e = 


V^VK 

f ■ 


(99) 


The energy density e of the ground state of the system is obtained by extremizing with respect to x, p, and 9. 
The pressure of the system is then obtained from the energy density 


P = p 


djs/p) 

dp 


( 100 ) 


The energy density e and pressure P are then used in the TOV equation to obtain the properties of neutron stars, 

dM 
dr 

dP _ G[e + P][l + 47rr3p] 
dr r{r — 2GM) 


= Anr^e 


( 101 ) 


In Fig. 48 we show the electron chemical potential as a function of nucleon density. In the case without kaons, 
the electron chemical potential continues to increase with nucleon density. When kaons are included, the electron 
chemical potential first increases up to about three times normal nuclear matter density po- Afterwards it starts to 
decrease because of the onset of the kaon condensation. Thus based on the kaon in-medium properties as constrained 
by heavy-ion data and nuclear equation of state from the model of Furnstahl, Tang and Serot, we find that the critical 
density for kaon condensation to be about 3po- 

In Fig. we show the proton fraction x and rotation angle 0 as a function of nucleon density. The solid and dotted 
curves are obtained with and without kaons. The proton fraction increases rapidly to become almost as large as the 
neutron fraction, after kaon condensation sets in. This reduces the asymmetry parameter, and hence the symmetry 
energy, of the system. The rotation angle, or the condensate amplitude, rises rapidly near the threshold, and then 
slowly increases to about 60°. 

In Figs. ^ and ^ we show the energy density and pressure of system as a function of nucleon density, for both 
cases with and with kaons. For densities lower than 0.08 fm“° we use the empirical equation of state as in Ref. 

The match between this and the equation of state of Furnstahl et al. produces the small dip near O.Spo- Both the 
energy density and pressure are seen to be lowered once kaon condensation sets in. The effect is particularly strong 
in the pressure. Since the properties of neutron stars also depend sensitively on the ratio P/e, we show in Fig. the 
pressure as a function of energy density. The softening of the equation of state by kaon condensation is appreciable. 

The results for neutron star mass as a function of central density Pcent are shown in Fig. |^. It is seen that, 
without kaons, the maximum neutron star mass in this model is about 2Mq. Similar conclusions, with neutron star 
mass in the range of 2.1-2.3M0, have been obtained in Ref. based on the nuclear equation of state from the 

Dirac-Brueckner-Hartree-Fock approach. When kaons are included, the maximum mass of neutron stars is about 
1 .5M0. The dependence of neutron star mass on radius is shown in Fig. It is seen that with kaon condensation, 
the radius of maximum mass neutron star increases. 


V. SUMMARY AND OUTLOOK 

In summary, we studied K~^ and K~ production in Ni-|-Ni collisions at 1-2 AGeV, based on the relativistic trans¬ 
port model including the strangeness degrees of freedom. We found that the recent experimental data from the KaoS 
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collaboration are consistent with the predictions of the chiral perturbation theory that the K'^ feels a weak repulsive 
potential and K~ feels a strong attractive potential in nuclear medium. Using the kaon in-medium properties con¬ 
strained by the heavy-ion data, we have studied neutron star properties with and without kaon condensation. The 
maximum mass of neutron stars is found to be about 2 .OM 0 based on conventional nuclear equations of state obtained 
from the effective Lagrangian of Furnstahl et ai. This can be reduced to about 1.5Mq, once kaon condensation is 
introduced. We have emphasized the growing interdependence between hadron physics, relativistic heavy-ion physics 
and the physics of compact stars in astrophysics. 

The recent experimental data from the COSY-11 collaboration on near threshold kaon and antikaon produc¬ 

tion cross section in proton-proton collisions provide much needed information for the parameterizations of elemen¬ 
tary cross sections. It is realized that both the Randrup-Ko parameterization for production cross section and 
Schiirmann-Zwermann parameterization for the K~ production cross section overestimate the data near threshold by 
a large factor. Use of the new parameterizations that incorporate the latest COSY-11 data reduces the contributions 
to and K~ yields in heavy-ion collisions from baryon-baryon interactions |^,101|. As a result, contributions 
from pion-baryon and pion-hyperon interactions, that were previously found to be less important than those from 
baryon-baryon interaction, turn out to be important, and in some cases even become dominant. 

Our knowledge on elementary cross sections has been improved over the last a few years, mainly because of the new 
experimental data near the threshold, and also because of the development in the theoretical study of the elementary 
cross sections. Nevertheless, to extract kaon medium effects from their yields in heavy-ion collisions is still a delicate 
task. On the other hand, the ratio, K^/K~, and especially the shapes of its transverse mass spectra and rapidity 
distribution, is less sensitive to the elementary cross sections, and should be able to p rovide m ore definite information 
about kaon in-medium properties. Also, the collective flow signals of and K~ 1102-104] in heavy-ion collisions 
are very useful probes of kaon medium effects. 
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FIG. 2. Kaon and antikaon potential in nuclear medium at two different momenta. The open circle in th e figur e 
potential expected from the impulse approximation using the kaon-nucleon scattering length in free space l48|J 
circle is the K~ potential extracted from the kaonic atom data ||l7|. 
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FIG. 3. Isospin-averaged cross sections a-rTN^YK- Solid and dotted lines are from Refs. and respectively. 
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FIG. 4. Kaon production cross sections ■ Open circles are experimental data, and the lines are parameterizations. 
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FIG. 5. Same as Fig. ^ for aTrN^SnK. 
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FIG. 6. Kaon production cross section iypp^pAK+ ■ The solid curve is the results of boson-exchanRe model of Ref. [^. The 
dashed and dotted curves are the parameterizations of Randrup and Ko and Gassing et al. [Q, respectively. The open 
circles are experimental data from Ref. while the solid circle is from Ref. [Q. 
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FIG. 7. Isospin-averaged kaon production cross sections (jbb^nyk ■ The solid lines are the results of boson-exchange model of 
Refs.^^^,^ . The dotted lines are the parameterizations of Randrup and Ko |Q . The dashed curves are the parameterizations 
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FIG. 8. Kaon production cross sections (Jnn^nKttK- The open circles are the experimental data from Ref. |^, and the 
lines are parameterizations. 
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FIG. 9. Same as Fig. ^ for onn^nt^ttk- 
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FIG. 10. Same as Fig. ^ for aNN^NA-nnK ■ 
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FIG. 11. Same as Fig. || for onn^nt^ttivk. 
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FIG. 12. Inclusive K'^ production cross section in pp collisions. The solid curve is the sum of pp —> NYK'^, NYttK^ , and 
NYmvK+. 
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FIG. 13. Antikaon production cross sections <^nN^NKK ■ The open circles are experimental data from Ref. and the lines 
are the parameterizations of Sibirtsev et al. [|^. 
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FIG. 14. Antikaon production cross section The open circles are experimental data from Ref. |^, and the 

lines are the parameterizations. 
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FIG. 15. Antikaon production cross section (^nn^nnkk- The open circles are experimental data from Ref. and the 
curves are our parameterizations. 
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FIG. 16. Isospin-averaged antikaon production cross section o-nn^nnkk ■ The solid line is from this work, and the dotted 
line is the parameterization of Zwermann and Schiirmann fiot. 
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FIG. 17. Same as Fig. E for Crjvjv-.iVArirKif • 
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FIG. 18. Same as Fig. ^ for o-ffN^NN-nT^KK- 
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FIG. 19. Inclusive K production cross section in pp collisions. The solid curve is the sum of pp —> ppK'^K , NNnKK'^ 
and NN'ktiKK~ . 








Pi.b (GeV/c) 

FIG. 21. Cross sections n^-jv^yTr- The open circles are experimental data from Ref. jM. The solid lines are the parame- 
terizations in this work, and the dotted lines are the parameterizations of Cugnon et al. [[f^. 
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FIG. 22. Cross sections a^y^KN- 
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FIG. 23. KN elastic, one-pion production, and two-pion production cross sections, 
from Ref. [^, and the curves are parameterizations. 


The symbols are experimental data 
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FIG. 24. K p total, elastic, and charge-exchange cross sections. The symbols are experimental data from Ref. p^ , 
curves are parameterizations. 
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FIG. 25. Same as Fig. for K n total and elastic cross sections. 
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FIG. 26. Angular distribution in tt p ^ AK° at = 1.742 GeV. The experimental data are taken from Ref. fisl ], while 
the solid curve is our parameterization. 
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FIG. 27. Kaon momentum spectra in pp collisions at = 2.88 and 2.99 GeV. The experimental data are taken from Ref. 
[jrit, while the the curves are parameterizations. 
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FIG. 28. Nucleon scalar and vector potentials from the effective chiral Lagrangian of Furnstahl, Tang and Serot [pl|. 
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FIG. 29. Proton transverse mass spectra in central Ni+Ni collisions. The open circles are experimental data from Ref. 
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FIG. 30. Same as Fig. 


for TT 
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FIG. 31. Proton rapidity distribution in central Ni+Ni collisions. The solid circles are experimental data from Ref. 
open circles are reflected with respect to mid-rapidity. 
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FIG. 32. Same as Fig. |3^, for tt . 
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FIG. 33. tranverse mass spectra and rapidity distribution in Ni+Ni collisions at 1.8 AGeV and 0 fm. The solid and 
dotted lines are the results with isotropic and anisotropic angular distributions, respectively. 
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FIG. 34. tranverse mass spectra and rapidit y d istribution in Ni+Ni collisions at 1.8 AGeV and 0 fm. The solid and 
dotted lines are the results using Eq. (M) and Eq. ([8^ as momentum spectra, respectively. 


57 






























104 


1000 ^ 


CO 

o 


CV2 

> 

cu 

O 


100 

CO 

a 

Xi 


b 

CO 

X! 

W 


10 ^ 


T 


T 


T 


1 —^ 




Ni + Ni, 0.8 AGeV, K+, 

.no kaon medium effects 

—with kaon medium effects 


C) i 


C 

%) 


-LC) 


5 


c:> 




4i 


1 


1 


1 


0 0.1 0.2 0.3 0.4 

EL" (GeV) 


0.5 0.6 


FIG. 35. kinetic energy spectra in Ni+Ni collisions at 0.8 AGeV. The solid and dotted histograms are the results with 
and without kaon medium effects, respectively. The open circles are experimental data from Ref. p^. 
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FIG. 36. Same as Fig. M, for Ni+Ni at 1.0 AGeV. 
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FIG. 37. Same as Fig. M, for Ni+Ni at 1.8 AGeV. 
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FIG. 38. Same as Fig. for K spectra in Ni+Ni at 1.8 AGeV. 
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FIG. 39. kinetic energy spectra of /K~ in Ni+Ni collisions at 1.0 AGeV {K^) and 1.8 AGeV {K~). The solid and dotted 
histograms are the results with and without kaon medium effects, respectively. The circles are the experimental data from the 
KaoS collaboration 


62 











































0.0001 


0 0.1 0.2 0.3 0.4 0.5 0.6 

m^-iTio (GeV) 


FIG. 40. K'^ transverse mass spectra in central Ni+Ni collisions at 1.0 AGeV. The solid and dotted histograms are the 
results with and without kaon medium effects, respectively. 
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FIG. 41. Same as Fig. ^ for central Ni+Ni collisions at 1.93 AGeV. 
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FIG. 42. Same as Fig. ^ for K in central Ni+Ni collisions at 1.8 AGeV. 


65 
















































dN/dy 



y 


FIG. 44. Same as Fig. ^ for central Ni+Ni collisions at 1.93 AGeV. 
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FIG. 45. Same as Fig. ^ for K in central Ni+Ni collisions at 1.8 AGeV. 
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FIG. 46. Rapidity distribution of jK in Ni+Ni collisions at 1.0 {K'^) and 1.8 {K ) AGeV. The solid and dotted curves 
are the results with and without kaon medium effects, respectively. 
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FIG. 47. Excitation function of jApart and jApart in central Ni+Ni collisions. The experimental data are from Ref. 
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FIG. 48. Electron chemical potential as a function of nucleon density. The solid and dotted curves are obtained with and 
without kaons. 
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FIG. 51. Pressure as a function of nucleon density. The solid and dotted curves are obtained with and without kaons 
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FIG. 52. Pressure as a function of energy density. The solid and dotted curves are obtained with and without kaons. 
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